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I.  INTRODUCTION 


An  initially  flat  circular  membrane  forming  one  end  of  a  cylindrical  vacuum  chamber,  and  clamped  along  the 
rim,  deforms  under  vacuum  to  a  surface  whose  shape  has  been  found  to  be  less  than  ideal  for  optical  applications. 
Deviations  from  an  ideal  paraboloid  or  sphere  are  greater  than  100A  at  optical  wavelengths  A,  even  for  membrane 
(paraxial)  /-numbers  as  high  as  //8.  The  work  reported  here  was  carried  out  in  support  of  an  in-house  laboratory 
effort  to  demonstrate  the  improvement  of  this  shape  by  introducing  a  controllable  strain,  independent  of  the  vacuum 
pressure  used  to  produce  the  membrane  curvature.  These  experiments  were  performed  on  a  Upilex  membrane  mounted 
on  two  concentric  cylinders  to  form  one  end  of  two  independently  controlled  vacuum  chambers.  When  a  vacuum  is 
created  in  the  interior  of  the  outer  annulus  (p0  <  patm ,  where  patm  is  laboratory  atmospheric  pressure),  with  no 
vacuum  on  the  inner  cylinder  (p  =  patm ),  the  membrane  is  free  to  move  outward  across  the  inner  ring,  creating  a 
strain  in  the  inner  membrane.  For  a  given  vacuum  pressure  p  <  patm  on  the  inner  cylinder,  this  “pre-strain”  favorably 
affects  the  final  shape  of  the  inner  membrane  surface  compared  to  the  shape  it  would  attain  with  no  pre-strain  (note 
that  when  p  <  patm,  the  annulus  pressure  may  change  to  some  new  value  p'0  due  to  migration  of  membrane  material 
back  across  the  inner  ring).  The  amount  of  pre-strain  can  be  controlled  by  the  annulus  pressure,  affording  a  degree  of 
control  over  the  final  shape  of  the  inner  membrane.  Details  of  this  “rim-controlled”  membrane  and  its  use  in  a  simple 
optical  system  were  presented  in  references  [6],  [7],  and  [9].  The  design  of  the  membrane  mounting  is  illustrated  in 
Figure  1. 
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FIG.  1.  In-House  Laboratory  Membrane  Mirror 
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The  goals  of  the  modelling  effort  described  in  this  report  were  to  provide  theoretical  verification  of  pre-strain 
control  of  a  membrane  surface,  and  to  investigate  effects  of  the  control  parameters  (inner  membrane  pressure  p  and 
pre-strain  e)  on  the  final  shape.  The  boundary  value  problem  corresponding  to  the  rim-controlled  membrane  mirror 
used  in  our  laboratory  experiments  is,  however,  rather  formidable,  and  to  the  author’s  knowledge  its  solution  has 
not  been  given  in  the  literature.  We  have  thus  opted  to  instead  examine  a  different  boundary  value  problem  whose 
solution  is  fairly  well-known,  and  which  corresponds  to  an  experimental  situation  similar  in  some  respects  to  the  one 
in  our  laboratory.  Results  obtained  from  this  solution  will  be  shown  to  be  in  qualitative  agreement  with  measurements 
made  on  the  laboratory  membrane,  and  can  thus  serve  as  a  guide  in  optical  designs  using  the  membrane  as  a  mirror. 
In  addition,  the  next  generation  of  experiments  at  this  laboratory  are  expected  to  involve  a  membrane  mounting  that 
more  nearly  resembles  the  one  described  by  the  boundary  value  problem  solved  here. 

As  a  quick  review  of  the  history  of  this  problem,  we  note  that  in  1915  Hencky  [5]  proposed  a  system  of  equations  for 
determining  the  equilibrium  configuration  of  an  initially  plane  circular  membrane  deformed  by  an  axially  symmetric 
constant  pressure  load.  They  are  essentially  von  Karman’s  nonlinear  plate  equations  [11]  specialized  to  zero  flexural 
rigidity  D  (for  a  clear  and  careful  discussion  of  the  von  Karman  equations,  see  [3]).  Hencky  derived  power  series 
solutions  for  the  stress,  strain,  and  displacement  vector  components  for  the  case  of  an  initially  unstressed  membrane. 
In  1956  Campbell  [1]  generalized  Hencky’s  problem  to  include  an  initial  tension  in  the  membrane.  In  Sections  II-V 
we  give  a  detailed  description  of  the  setting  for  the  Hencky-Campbell  problem,  and  show  that  it  is  equivalent  to  the 
problem  of  a  membrane  subjected  to  an  initial  purely  radial  displacement  of  its  edge,  then  clamped  at  the  original 
radius  before  being  deformed  by  a  uniform  pressure.  This  boundary  value  problem  is  illustrated  in  Figure  2,  where 
the  differences  between  the  laboratory  membrane  mounting  and  the  one  modeled  by  the  Hencky-Campbell  theory 
are  apparent,  i.e.,  the  Hencky-Campbell  membrane  is  clamped  at  the  inner  ring  after  pre-straining  it  by  a  radial 
displacement  of  its  edge,  and  is  then  pressurized,  while  the  laboratory  membrane  can  slip  back  over  the  inner  ring 
when  it  is  pressurized,  after  pre-straining  it  by  the  annulus  pressure. 
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FIG.  2.  Hencky-Campbell  Membrane  Model 
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Sections  VI  and  VII  are  concerned  with  the  elasto-mechanical  properties  of  the  membrane  derived  from  the  Hencky- 
Campbell  theory,  while  Section  VIII  begins  a  discussion  of  relations  between  these  properties  and  some  of  the  optical 
parameters  of  the  problem. 


II.  REFERENCE  CONFIGURATION 

We  consider  a  body,  henceforth  referred  to  as  a  membrane,  whose  reference  configuration  is  a  thin  circular  disk  of 
radius  a  and  thickness  h  «  a.  An  arbitrary  point  P  of  the  reference  configuration  can  be  specified  by  either  its 
Cartesian  coordinates  {X,  Y,  Z},  or  its  cylindrical  coordinates  (fi,  Q,Z).  The  middle  plane  of  the  membrane  is  the 
plane  Z  -  0,  and  we  take  the  origin  of  coordinates  to  be  the  center  0  of  the  circle  of  radius  a  in  this  plane.  We 
assume  given  a  fixed  orthonormal  Cartesian  basis  {i,j,fc}  at  0 .  The  position  vector  of  P  with  respect  to  O  is  given 
by 

X  =  Xi  4  Yj  4  Zk.  (2.1) 

The  cylindrical  and  Cartesian  coordinates  of  P  are  related  by  the  following  transformation  equations: 

X  —  R  cos0,  Y  —  R  sinG,  Z  =  Z,  (2.2) 

hence  the  position  vector  can  be  written  in  terms  of  cylindrical  coordinates  as 

X  =  RcosQi  4  RsinQj  4  Zk.  (2*3) 

We  introduce  orthonormal  basis  vectors  {Er,  25©,  Ez}  associated  with  the  cylindrical  coordinates,  and  defined  in 
terms  of  the  fixed  basis  vectors  by 

Er(Q)  =  cos  ©  i  4  sin  ©  j,  15©  (0)  =  -sin0i  +  cos©  j,  Ez  =  k,  (2*4) 

in  terms  of  which  the  position  vector  can  be  written  as 

X{R,e,Z)  =  REr(Q)  4  ZEz.  (2-5) 


III.  ASSUMPTIONS  ON  THE  DEFORMATION 

Under  a  uniform  pressure  p,  a  point  P  of  the  reference  configuration  is  deformed  to  a  new  point  P  .  The  set  of  all 
such  points  defines  the  deformed  configuration.  The  point  P'  of  the  deformed  configuration  has  Cartesian  coordinates 
z}  and  cylindrical  coordinates  (r,  0,z).  The  position  vector  of  Pf  with  respect  to  the  origin  O  at  the  center  of 
the  membrane  is  thus  given  in  the  fixed  Cartesian  basis  by 

x  —  xi  yj  +  zk  =  —  r  cos#«  4-  r  sin 9  j  4  zk,  (3*1) 

where  the  cylindrical  and  Cartesian  coordinates  are  related  by 

x  =  r  cos#,  y  —  r  sin0,  z  =  z.  (3.2) 

The  orthonormal  basis  vectors  {er,  ee,ez}  associated  with  the  cylindrical  coordinates  on  the  deformed  configuration 
are  defined  by 

er(0)  =  cos 6i  +  sin#  j,  ee{6)  =  —sin Oi  4  cos 9 j,  ez  =  k  =  Ez ,  (3*3) 

in  terms  of  which  the  position  vector  can  be  written  as 

x  =  rer4^e:.  (3*4) 

The  displacement  vector  u  relates  the  position  vector  X  of  P  in  the  reference  configuration  to  the  position  vector 
x  of  the  point  Pf  to  which  P  is  displaced  by  the  deformation,  i.e.,  u  =  x  -  X,  hence 

*  =  X  4  w-  (3-5) 
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We  write  the  displacement  vector  in  the  cylindrical  basis  {Er,E©,  Ez}  of  the  reference  configuration  as 

u  =  UrEr  4  U©E©  4  UzEz,  (3.6) 

where  the  orthonormal  cylindrical  basis  vectors  are  related  to  the  fixed  Cartesian  basis  vectors  by  (2.4).  Thus,  (3.5) 
can  be  written  as 

xi  4  yj  4  zk=z  Xi  +  Y  j  4  Z  k  4  Ur  Er  4  U ©  E©  +  Uz  Ez 

=  Xi  4  Y  j  +  Zk  4  Ur  ( cos  0  i  4  sin  0  j  )  +  U©  (  —  sin  0  i  4  cos  Qj)  4  Uz  k, 

—  ( X  4  Ur  cos  0  —  U©  sin  0 )  i  4  ( Y  4  Ur  sin  ©  4  U©  cos  0 )  j  4  ( Z  4  Uz  )  (3.7) 

From  (3.7)  we  obtain  the  following  relations  between  the  Cartesian  coordinates  of  the  two  position  vectors,  and  the 
cylindrical  components  of  the  displacement  vector: 

x  =  X  4  Ur  cos©  -  U©  sinG,  (3.8a) 

y  =  y  4  Ur  sin©  4  U©  cos©,  (3.8b) 


z  =  Z  4  Uz-  (3.8c) 

Substituting  from  (2.2)  and  (3.2)  into  the  first  two  equations  of  (3.8),  we  obtain 

r  cos#  =  (R  4  Ur)  cos  ©  —  U©  sin©,  (3.9a) 

r  sin#  =  (R  4  Ur  )  sin  0  4  U©  cos©.  (3.9b) 

These  equations  can  be  solved  for  the  deformed  configuration  cylindrical  coordinates  in  terms  of  the  reference 
configuration  cylindrical  coordinates  and  displacement  vector  components  as  follows:  first,  multiply  the  first  equation 
by  cos©  and  the  second  by  sin©,  and  add  the  results  to  get  (using  the  trigonometric  identity  for  the  cosine  of  the 
difference  of  two  angles)  rcos  (0  -  6)  =  R  4  Ur;  second,  multiply  the  first  equation  by  -  sin©  and  the  second  by 
cos©,  and  add  to  get  rsin  (0  —  6)  =U©.  From  these  two  results  one  easily  obtains 

r  =  y/(R+  t/fl)2  +  Ul,  ,  (3-lOa) 

=  ,310b) 

We  assume  that  the  deformation  of  the  membrane  from  its  reference  configuration  to  its  deformed  configuration  is 
described  mathematically  by  one-to-one  invertible  mappings  between  the  coordinates  of  points  in  the  two  configura¬ 
tions,  e.g.,  between  their  cyindrical  coordinates: 

r  =  fr(R,e,Z),  0  =  fe(R,  0, Z),  z  =  fz(R,e,Z).  (3.11) 

We  immediately  restrict  consideration  to  deformations  in  which  the  cylindrical  coordinate  ©  of  any  point  in  the 
reference  configuration  remains  the  same  after  the  deformation,  i.e.,  we  assume  that  the  second  equation  of  (3.11)  is 
the  identity  mapping 

6  =  0.  (3.12) 

Furthermore,  we  restrict  consideration  to  deformations  that  are  independent  of  this  angular  coordinate.  Equations 
(3.11)  thus  reduce  to  a  pair  of  mappings  depending  only  on  the  other  two  coordinates: 

r  =  fr(R,Z),  z  =  fz(R,  Z).  (3.13) 

Under  these  assumptions  of  axisymmetric  deformation,  equations  (3.12)  and  (3.10b)  imply  that 

Ue  =  0,  (3.14) 
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hence  the  relation  (3.10a)  between  radial  coordinates  in  the  two  configurations  simplifies  considerably  to 

r  =  fr{R,  Z)  =  R+  Ur{R,  Z).  (3.15) 

The  relation  (3.8c)  between  the  vertical  coordinates  is  unchanged: 

z  =  fz(R,Z)  =  Z  +  Uz(R,Z),  (3.16) 

and  the  displacement  vector  reduces  to 

u(R,e,Z)  =  Ur{R,Z)Er{Q)  +  Uz{R,Z)Ez.  (3.17) 

Assumption  (3.12)  implies  that  the  orthonormal  basis  vectors  {eR,e@,ez}  on  the  deformed  configuration  are  the 
same  as  those  on  the  reference  configuration,  i.e., 

er(9)  -  er(0)  =  Efi(0),  es(9)  -  e<>(@)  =  E@(Q),  ez  =  Ez  =  k.  (3.18) 


IV.  HENCKY’S  EQUATIONS 

Hencky’s  equations  describe  the  equilibrium  state  of  the  middle  surface  of  the  deformed  membrane,  i.e.,  the  image 
under  the  deformation  of  the  middle  plane  Z  =  0.  Setting  Z  =  0  in  (3.15)  and  (3.16),  we  write  the  equations  describing 
the  middle  surface  as 

r  =  F(R)  =  R  +  U(R),  where  F{R)  =  /r(-R,  0)  and  U(R)  =  C/r(.R,0),  (4.1) 

and 

z  =  W{R),  where  W(R)  =  f,{R,  0)  =  UZ(R,  0).  (4.2) 

The  deformation  of  the  middle  plane  into  the  middle  surface  is  illustrated  in  Figure  3,  below. 


We  consider  a  circular  membrane  made  of  a  material  characterized  by  given  values  of  Young’s  modulus  E  and 
Poisson’s  ratio  v.  We  assume  that  under  a  uniform  pressure  p  the  membrane  is  deformed  to  a  new  equilibrium 
configuration,  illustrated  in  Figure  3,  that  is  a  solution  of  Hencky’s  equations  [5].  This  system  of  equations  consists 
first  of  the  following  relations  between  the  radial  and  circumferential  strain  tensor  components  and  the  displacement 
vector  components: 


_  dU_ 
~  dR 


1  rdW\ 2 

2  \dR)  ’ 


(4.3) 
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(4.4) 


U 

£Q  =  R' 

where  all  variables  are  assumed  to  depend  only  upon  the  radial  coordinate  R  of  the  reference  configuration.  The 
strain  components  are  assumed  to  be  linearly  related  to  the  second  Piola-Kirchhoff  stress  components  (see  [3],  pp. 
465-466)  by  the  following  uniform,  homogeneous,  isotropic  form  of  Hooke’s  Law: 


Sr  -  v  S®  —  EeR , 

(4.5) 

5©  —  V  Sr  —  E  6© . 

(4.6) 

Finally,  the  stress  components  must  satisfy  the  force  equilibrium  equations 

50  =  SR  +  R ^  =  Jn(RSR), 

(4.7) 

dW  _  _pR 
Sr  dR  ”  2h  ’ 


(4.8) 


representing  equilibrium  in  the  radial  and  lateral  directions,  respectively.  We  remark  that  Fichter  [2]  has  recently 
suggested  that  (4.8)  represents  equilibrium  in  the  lateral  direction  only  for  a  laterally  applied  pressure  load,  and  should 
be  altered  to  include  an  additional  term  if  the  model  is  to  represent  a  true  pressure  load  (which  would  necessarily 
be  normal  to  the  deformed  surface).  Regardless  of  the  merits  of  his  suggestion,  which  at  first  sight  appear  to  be 
well  motivated  physically,  we  adopt  in  this  report  the  more  familiar  model  using  the  equilibrium  equation  (4.8).  We 
believe  that  Fichter’s  claim,  if  true,  should  be  derived  in  a  way  that  carefully  distinguishes  between  variables  defined 
on  the  reference  and  deformed  configurations.  Such  a  derivation  has,  to  our  knowledge,  not  been  given. 

Before  going  on  to  consider  a  particularly  simple  solution  of  these  equations,  we  give  details  of  the  derivations  of 
two  useful  differential  equations  that  follow  from  (4.3)-(4.8).  Differentiating  (4.4),  and  using  (4.3)  and  (4.4)  in  the 
result,  we  obtain 


which  can  be  written  as 


de©  1  dU 
dR  R  dR 


U 

R? 


r  i 

fdW\*-\ 

b- 21 

Iir)  J 

=  tR 


C© 


1  /dW\ 2 

2  \dRJ  ' 


£© 

R' 


Multiplying  this  equation  thru  by  E  and  using  (4.5)  and  (4.6)  to  replace  the  strains  in  terms  of  the  stresses  yields 


d  E  fdW\ 2 

RdR{S@  -  vS*)  =  -  uRe  ~  S&  +  vS*  ~  J  VdRJ  ' 

On  the  right-hand  side  of  the  last  equation,  we  use  (4.7)  to  replace  Sr  =  5©  —  R(dSR/dR),  obtaining  after  some 
simplification: 


RjR{Se  +  Sr)  +  f  (S)  “  °’  (4'9) 

which  is  Hencky’s  equation  (4),  and  Campbell’s  equation  (7).  This  can  be  written  in  terms  of  the  pressure  by  using 
(4.8)  to  replace  dW/dR ,  yielding 


tr^+s^+ !£=»• 


It  is  convenient  to  introduce  at  this  point  a  dimensionless  independent  variable  p  defined  by 


P  = 


R 
a  ’ 


(4.10) 


(4.11) 


as  well  as  dimensionless  displacement  vector  components  U  and  W,  defined  by 
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(4.12) 


U{p)  =  \u{ap),  W{p)  =  \w{ap ), 
and  dimensionless  stress  components  Sr  and  5© ,  defined  by 

Sr{p)  =  “  SR(ap ),  Sq{p)  =  ~  Se(ap)-  (4-13) 

We  have  departed  somewhat  from  the  Hencky- Campbell  approach,  as  they  define  dimensionless  stress  components 
in  terms  of  a  factor  k  that  is  proportional  to  p.  With  our  definitions,  equations  (4.3)-(4.8),  (4.9),  and  (4.10)  can  be 
written  in  dimensionless  form  as 


dU  1  (dW\2 
€R  ~  dp  +  2  \  dp  )  ’ 

(4.14) 

U 

e©  —  — , 

P 

(4.15) 

Sr  —  V  Sq  —  €r, 

(4.16) 

Sq  —  v  Sr  =  e©, 

(4.17) 

5©  =  Sr  +  p  ^  =  —  (pSn), 

(4.18) 

5  dW  q 

Sr1^  =  -2P’ 

(4.19) 

and 


respectively,  where 


+  «•>  +  *  (fr-' 

-f(Se  +  SR)  +  \q2  =  0, 

p  dp  o 


is  a  dimensionless  constant  proportional  to  the  external  pressure  p. 


(4.20) 

(4.21) 

(4.22) 


V.  EQUILIBRIUM  CONFIGURATION  OF  AN  UNPRESSURIZED  CIRCULAR  MEMBRANE 

We  consider  here  the  solution  of  equations  (4.14)-(4.19)  when  there  is  no  externally  applied  pressure,  which  must 
include  as  one  case  the  initially  tensioned  state  assumed  in  Campbell’s  paper.  With  the  exception  of  the  radial 
coordinate  on  the  deformed  configuration,  all  variables  in  this  special  situation  will  carry  a  0-subscript,  hence  the 
deformation  equations  corresponding  to  (4.1)  and  (4.2)  will  be  written  as 

R  =  R0  -(-  Uq(Ro)  and  z0  =  Wo(J?o),  0>-l) 

respectively.  We  also  introduce,  similarly  to  (4.11),  a  dimensionless  radial  coordinate  p0  defined  by 
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o 

III 

(5.2) 

Thus,  #o  (or  po)  labels  a  point  of  the  undeformed  reference  configuration,  and  R  (or  p)  labels  the  point  to  which  it 
may  be  deformed  in  the  absence  of  external  pressure. 

It  is  convenient  to  again  introduce  dimensionless  displacement  vector  and  stress  tensor  components  defined  by 

Uo{po)  =  -Uo{apQ),  Wo(po)  =  -  W{apv), 

a  cl 

(5.3) 

and 

Sro(po)  =  -gr  SRo{ap0),  S@0{po)  =  -j=r  Seo{apo)- 

(5.4) 

From  (5.3),  (4.11),  (5.2)  and  (5.1),  we  obtain  a  dimensionless  form  of  the  equation  of  radial  deformation: 

p  =  po  -f  Uo(po)* 

(5.5) 

The  system  of  equations  (4.14)-(4.19)  with  p  =  0  (hence  q  =  0  )  can  be  written 
variables  as 

in  terms  of  these  dimensionless 

_  dU0  1  (dWo\> 

€Ra)  dpo  2  1  dpo  ) 

(5.6) 

Uo 

e©o  =  — , 

Po 

(5.7) 

Sro  —  J'S'eo  =  €rq, 

(5.8) 

Se  o  —  v  Sro  =  e©o, 

(5.9) 

Seo  —  Sro  +  po  =  ,  ( PoSro ), 

dpo  dpo 

(5.10) 

Sm  f  =  0, 

apo 

(5.11) 

It  follows  immediately  from  (5.11)  that  either  Srq  is  zero,  or  Wo(po)  is  a  constant,  or  both  of  these  conditions  hold. 
It  is  easy  to  see  that  if  Srq  =  0,  then  all  dependent  variables  except  Wo  vanish,  and  Wq  is  again  relegated  to  some 
constant  value  by  equation  (5.6). 

The  remaining  possibility,  then,  is  that  Wo(po)  is  constant,  but  Sro  ^  0.  Under  these  conditions,  the  lateral 
equilibrium  equation  (5.11)  is  identically  satisfied,  and  equation  (5.6)  for  the  radial  strain  reduces  to 


_  dU0 
dp0 


(5.12) 


hence  the  system  we  must  solve  consists  of  equations  (5.7),  (5.8),  (5.9),  (5.10),  and  (5.12),  subject  to  the  requirement 
that  all  functions  be  regular  (i.e.,  be  defined  for  all  values  of  po).  We  begin  by  using  (5.10)  to  replace  S© o  in  the 
left-hand  sides  of  (5.8)  and  (5.9),  and  using  (5.7)  and  (5.12)  in  the  right-hand  sides  of  the  same  two  equations,  to 
write  them  as 


( 1  “  v )  Srq  —  v  po 


dSRo 

dpo 


dUp 

dpo  ’ 


and 
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( 1  —  v )  Srq  +  po 


dSno  Uq 


dpo  po 

Subtracting  the  first  of  these  two  equations  from  the  second,  we  obtain 

dSfio  /  Uq  dUo 


(1  +  v)po 


/  uo  \ 

dpo  '  po  dpo  ' 


which  can  be  rewritten  as 


dSm 

dpo 


_ _ 1_  ,  ]_d%  _  Uo  v  _ 1 _ <L(U±\ 

~  1  +  V  ^  Po  dpo  Po  ~  1  +  V  dpo  V  Po  )  ‘ 


This  integrates  immediately  to  obtain 


V  1  +  v  )  po  ’ 


Srq  =  ci 

where  C\  is  an  arbitrary  integration  constant.  Substitution  of  (5.13)  and  (4.4)  in  (4.6)  yields 

-  /  1  \  Uq 

Finally,  substitution  of  (5.13),  (5.14)  and  (5.12)  in  (5.8)  yields  a  differential  equation  for  Uo(po)- 


Cl 


-(rb) 


Uo 

+  V  )  PQ 


Cl 


1  \  Uq  1  dUo_ 

dpo  ’ 


,  f  Uq  1  _ 

V  \  1  4*  v  /  Po  -I 


which  can  be  rearranged  to  write  it  as 


which  easily  integrates  to  obtain 


dpo 


{po  Uq  )  =  Ci  ( 1  -  i/2)  po, 


(5.13) 


(5.14) 


Tr  /  X  _  1  -  *  ,  C2 

Uq(pq)  -  ci — - — po  +  —  ■ 
l  Po 

In  this  expression,  the  arbitrary  integration  constant  C2  must  be  set  equal  to  zero  in  order  for  Uq  to  be  regular  at 
pQ  =  0,  hence  the  general  solution  for  Uq  reduces  to 


^  1  —  i/2 

tM/>o)  —  ci — 2 — ^°* 

Substitution  of  this  result  in  (5.6)  and  (5.7),  recalling  that  W0  is  constant,  yields 


tRO  —  C0Q  —  Ci 


1  -  i/2 


(5.15) 


(5.16) 


i.e.,  the  two  strain  components  are  equal  to  the  same  constant  value,  which  we  denote  by  e.  From  (5.8)  and  (5.9)  it 
then  follows  that  the  two  stress  components  are  also  equal  and  constant  throughout  the  membrane: 


Srq  =  5©0  =  6  =  S. 


(5.17) 


In  terms  of  the  constant  strain  e,  the  radial  displacement  (5.15)  is  given  by  the  linear  relation 

Uo{po)  =  cpQ  =>  Uq(1)  =  e.  (5.18) 

Thus,  Campbell’s  configuration  of  initial  tension  a  (which  corresponds  to  5  =  E  S  in  the  preceding  development)  is 
equivalent  to  a  configuration  in  which  each  point  has  suffered  a  linear  radial  displacement  from  its  unstrained  reference 
configuration.  The  radial  deformation  (5.5)  is  found  from  (5.18)  and  (5.16)  to  have  the  simple  form 

P  =  po{  1  +  0,  (5.19) 

while  the  lateral  deformation  is  taken  to  be 

zq  =  0,  (5.20) 

assuming  the  constant  lateral  displacement  to  be  zero,  i.e.,  Wo(po)  —  0.  Since  Wq(pq)  is  constant,  this  will  be  true  if, 
for  example,  the  membrane  remains  in  contact  with  the  hoop  edge  at  Rq  =  a  during  the  radial  displacement. 
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VI.  CAMPBELL’S  PROBLEM 


We  consider  now  a  membrane  which  has  suffered  a  purely  radial  displacement  of  its  edge.  As  shown  in  the  previous 
Section,  this  introduces  an  initial  uniform  dimensionless  stress  S  and  corresponding  strain  e,  related  by  (5.17).  Body 
points  of  the  middle  plane  of  this  pre-strained  membrane  are  labeled  by  the  radial  coordinate  R,  as  in  the  previous 
Section.  The  membrane  is  subsequently  clamped  along  the  circle  R  =  a,  and  an  external  uniform  pressure  is  applied 
to  points  of  the  disk  R  <  a.  In  terms  of  dimensionless  variables  defined  at  points  of  this  pressurized  part  of  the 
pre-strained  membrane  labeled  by  the  dimensionless  radial  variable  p  =  Rj a,  the  Hencky-Campbell  equations  can  be 
written  as 


dU  1  pdW\2 

CR  -  €  +  Tp  +  2  \~d^j  ’ 

(6.1) 

u 

c©  —  e  H - , 

P 

(6.2) 

Sr  —  vSq  =  £r, 

(6.3) 

1 

8° 

II 

© 

(6.4) 

(jl  d  ^ 

5©  =  Sr  +  P—(Sr)  =  —  ( P§r ), 

(6.5) 

$  dW  q 

SR  dp  2P’ 

(6.6) 

+  +  )!  =  °' 

(6.7) 

and 

ls^rSe  +  SR)  +  \q^0, 

(6.8) 

We  remark  that  the  stress  and  strain  components  in  these  equations  are  the  total  stresses  and  strains,  not  the 
incremental  stresses  and  strains  used  by  Campbell.  The  non-zero  pre-strain  is  accounted  for  in  equations  (6.1)  and 
(6.2),  which  are  modifications  of  Hencky’s  equations  (4.14)  and  (4.15)  needed  to  assure  that  the  jbotal  stress  and 
strain  in  the  membrane  reduce  to  the  pre-stress  and  pre-strain  when  the  displacement  components  U  and  W  vanish 
(corresponding  to  a  reduction  of  the  external  pressure  to  zero).  Since  the  membrane  is  clamped  along  the  circle  R  =  a, 
corresponding  to  p  =  1,  the  appropriate  boundary  conditions  to  be  applied  are 

U{  1)  =  0,  and  W(l)  =  0. 

(6.9) 

Similarly  to  Hencky  and  Campbell,  we  begin  by  assuming  an  even  power  series 
stress  component  Sr  of  the  form 

solution  for  the  dimensionless  total 

Sr  =  bo  ( 1  +  &2 nP2n  )i 

n=l 

(6.10) 

where  the  coefficients  ft  —  1,2,3, .. .,  are  dimensionless  constants.  Substitution  of  this  expression  in  (6.5)  yields 

5©  =  —  (p Sr)  =  bo  [  1  -h  ^2  +  l)&2np2n  ] • 

"  n  =  l 

(6.11) 

10 


From  the  last  two  equations  we  obtain 

Sr  +  S©  =  2  6q  [1  +  ^  (n  -f  l)b2 nP2n  ], 

hence 


n  =  1 


—  —  (  Sr  +  5©  )  =  4  6o  n(n  4-  l)&2np 


,2n  —  2 


P  dp 


n=l 


This  is  substituted,  together  with  (6.10),  into  (6.8)  to  obtain 

46§  (i  +  £w2")2  •£«("+ i)w2n~2  +  =  °> 

n= 1  n  =  l 

which  we  rewrite  as 

( 1  +  £  )2  ■  2  +  8!  =  0,  (6.12) 
n= 1  n=l 

where  we  have  introduced  a  new  dimensionless  constant  B  defined  by 

»-i£  *  -=i(sr  <6i3> 

°o 

The  expanded  form  of  (6.12)  is 

( 1  +  b2p2  +  &4p4  +  bep6  -f  ...  )2  (62  +  364/?2  4-  6£>6P4  4-  10&8P6  +  •  •  • )  +  52  =  0.  (6-14) 

This  is  used  to  obtain  all  coefficients  b2ni  n  >  1,  in  terms  of  B,  by  computing  the  products  of  the  three  infinite  series, 
and  equating  to  zero  each  coefficient  of  different  powers  of  p.  We  have  made  use  of  the  Mathematical  software  system 
to  determine  that  each  of  these  coefficients  has  the  form 


b2n  =  “A  nS2", 

where  the  purely  numerical  coefficients  /32n  for  n  =  1  to  n  =  11  are  given  by 


(6.15) 


A  =  li  A  =  A  —  jgi 


A  -  18i  Ao  27’ 


A  2  = 


1205 


_  219241 

Pi  4  —  /.OErtil  5  Pi  6 


6634069 


567  ’  ™  "  63504  ’  Ki0  1143072’ 

a  51523763  a  998796305  0  _  118156790413 

P 18  ClilOOOyl  5  /^20  r /»conA/?  >1  *  '22 


5143824  ’  56582064  ’  ^  3734416224  '  ^6'16^ 

Now,  substituting  (6.2)  in  the  right-hand  side  of  (6.4),  and  evaluating  the  result  at  the  edge  boundary  p  =  1,  yields 

60  ( 1  +  362  +  hb4  +  7&6  +  968  +  . . .)  -  vb0  ( 1  +  b2  +  b4  +  b6  +  b8  +  .  - .)  =  €  +  U(l)  =  e, 

where  we  applied  the  first  boundary  condition  (6.9).  This  can  be  rewritten,  after  collecting  like  terms  on  the  left-hand 
side,  as 


[(1  -  v)  +  (3  -  v)h  +  (5  -  v)b4  +  (7  -  u)be  +  (9  -  v)h  +  ••■]  -  ^  =  °- 


(6.17) 


Since  all  the  coefficients  have  been  expressed  in  terms  of  B,  or  equivalently,  in  terms  of  6o,  the  left-hand  side  of 
(6.17)  defines  a  function  whose  independent  variable  is  &o-  According  to  (6.17),  the  remaining  constant  b0  must  be 


*  Mathematica  is  a  registered  trademark  of  Wolfram  Research,  Inc. 
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determined  as  a  zero  of  this  function.  This  can  be  done  by  standard  numerical  methods,  e.g.,  the  bisection  method 
given  in  [8]  is  quite  satisfactory  for  our  purposes. 

Once  the  coefficient  bo  has  been  determined,  the  dimensionless  stress  components  are  given  by  the  even-order  power 
series  (6.10)  and  (6.11),  while  U(p)  is  given,  according  to  (6.2)  and  (6.4),  by  the  odd-order  series 

U(p)  =  p{Sq  —  v  Sr)  —  pe  =  bo  p[(  l  —  v)  +  ( 2  n  +  1  —  v)  b2nP2n]  —  pc.  (6.18) 

n  =  1 

Similarly,  the  two  strain  components  are  obtained  by  substituting  the  appropriate  series  into  (6.3)  and  (6.4),  to  obtain 

cr  =  Sr  —  v  Se  =  &o{(l  —  +  53  [1  —  J'' ( 2n -{- 1 )  ]  &2  nP2n)i  (6.19) 

71  =  1 

and 

e©  —  Sq  —  v  Sr  =  &o  {  ( 1  ~  v )  4"  5  ]  ( 2n  4-  1  —  ^ )  ^2n^2n  )•  (6.20) 

n  = 1 


The  remaining  dimensionless  variable  W{p)  can  be  determined  by  assuming  an  even-order  power  series  of  the  form 


W(p)  =  ^  w2„  p2n, 


n— 0 


the  derivative  of  which  is 


dW 

dp 


=  2  n  w2n  p2n  1 


n= 0 


Substitution  of  (6.22)  and  (6.10)  in  (6.6)  yields 


260  ( i  +  53  b^p2n )  53 n  ™2n  p2n  1 4"  | 9  =  °’ 


n  =  l 


n=0 


(6.21) 


(6.22) 


or 


( 1  +  62n/»2n )  ^2  nc2np2n  1  +  p  =  0. 


n=l 


n=0 


(6.23) 


(6.24) 


In  the  last  line  we  have  defined  new  coefficients  c2n  such  that  u>2n  =  C  c2n,  n  >  1 ,  where 

C  = 

4  60 

is  another  dimensionless  constant.  Expanding  the  product  of  the  two  infinite  series  in  (6.23),  replacing  each  &2n  in 
terms  of  B,  and  then  equating  coefficients  of  like  powers  of  p  to  zero  we  find,  again  making  extensive  use  of  the 
Mathematica  system,  that 

c2„  =  -72  nB2n~\  (6-25) 

where  the  purely  numerical  coefficients  j2n  for  n  =  1  to  n  =  10  are  given  by 

1  5  55  7 

72  =  1,  74  =  2  >  76  =  g,  78  =  710  =  6’ 


712  108’  714  5292’ 


718  = 


205 

108’ 
741805 

10287648’ 


17051  2864485 

716  = 


508032 


720  = 


16659221 
5143824  ’ 


(6.26) 


The  remaining  coefficient  w0  is  determined  by  the  second  boundary  condition  of  (6.9),  viz.,  W(l)  =  0,  which  yields 
from  (6.21): 

(6.27) 


0  =  W0  +  C  ^2  C2n  =>  Wp  —  —C 


C2n- 


n=l 


n  —  1 


Power  series  expressions  for  each  of  the  dimensionless  variables  are  now  completely  determined.  Each  of  these  must 
be  truncated  after  a  finite  number  of  terms,  determined  by  the  error  acceptable  by  the  user. 
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VII.  THE  EQUATION  OF  THE  “TRUE”  MIDDLE  SURFACE 

The  solutions  U(R)  and  W(R)  can  be  substituted  in  equations  (4.1)  and  (4.2)  to  obtain  the  deformation  mappings: 

r  =  F(R)  =  R  +  U(R),  and  z  =  W{R).  (7.1) 

They  can  be  interpreted  as  parametric  equations  of  the  middle  surface,  with  parameter  R.  The  second  equation  of  (7.1) 
defines  the  surface  of  lateral  displacements  of  the  body  points  with  radial  coordinates  R.  The  graph  of  this  equation 
is  frequently  identified  with  the  middle  surface,  but  it  does  not  define  the  true  middle  surface  (this  observation  has 
been  emphasized  in  a  recent  AIAA  Proceedings  paper  [4]  ).  Points  on  the  middle  surface  are  coordinatized  by  (■ r,z ), 
not  (R,z),  hence  the  equation  of  the  middle  surface  must  relate  the  lateral  displacement  z  to  the  radial  coordinate 
r  of  the  point  in  the  deformed  configuration  that  was  originally  at  R.  Thus,  it  must  take  into  account  the  radial 
displacement  U(R)  as  well.  As  can  be  seen  in  Figure  3,  the  true  middle  surface  should  be  graphed  by  plotting  the 
lateral  displacement  z  =  W  (R)  versus  the  point  of  the  middle  plane  with  coordinate  r  =  R  +  U(R)  {not  against  the 
point  with  coordinate  R  ). 

For  some  purposes  it  is  convenient  to  actually  eliminate  the  parameter  R  between  the  two  relations  in  (7.1).  We 
first  introduce  on  the  deformed  configuration  a  dimensionless  radial  variable  defined  by 


t  = 


(7.2) 


so  that,  recalling  the  definition  p  =  R/a ,  hence  R  =  ap,  we  can  evaluate  the  two  relations  in  (7.1)  at  R  —  ap  to  obtain 

£  =  -F(R)  =  -F(ap)  =  F(p)  =  p  +  U(p),  and  z  =  W(R)  =  W(ap )  =  aW{p),  (7.3) 

a  a 

where  we  recall  that  U (ap)  =  a  U(p).  One  can,  in  principle,  solve  the  first  of  these  relations  to  obtain  p  as  a  function 
of  £,  i.e., 


P  =  F-'it), 

which  can  then  be  substituted  into  the  second  relation  to  obtain 

z  =  W(R)  =  aW(F_1(£))  =  au;(£)  =  aw(r/a)  =  w(r). 


(7.4) 


(7.5) 


The  final  equality  is,  formally,  the  equation  of  the  true  middle  surface. 

Of  course,  it  is  difficult,  if  not  impossible,  to  determine  the  exact  relation  (7.4).  However,  such  a  relation  can  be 
found  to  any  desired  accuracy,  as  follows.  We  substitute  the  infinite  series  (6.18)  for  U(p)  into  the  first  relation  in 
(7.3),  and  write  the  result  as 

£  =  F(p)=  [(1  -  e)  +  60(1  -  v)]p  +  6062(3  -  v)  pz  +  b0b4(b  -  v)  p*  +  ..., 

=  d\  p  +  d$  p3  +  (is  p5  +  . . .  . 

We  then  express  the  variable  p  as  an  infinite  series  in  viz., 

p  =  f~HO  =  fit  +  fst3  +  fst5  +  •••  , 


(7.6) 


(7.7) 


where  the  f2n+i  are  determined  by  substituting  (7.7)  for  p  into  (7.6),  and  equating  like  powers  of  £  in  the  result.  For 
example,  we  find  in  this  way 


h  =  f3~ 


<h_ 

dV 


h 


3^  — 

a?  ’ 


h  — 


Sdid^ds  —  12d§  —  d^di 

dp  ’ 


(7.8) 


again  making  use  of  the  Mathematica  system.  In  practice,  the  coefficients  of  the  highest  order  term retained  in  the 
power  series  (7.6)  and  (7.7)  are  chosen  to  satisfy  the  boundary  condition  on  the  circle  p  =  1,  viz.,  U(l)  =  0,  which 
in  turn  requires  that  £  =  1  when  p  =  1.  For  example,  we  find  that  for  pre-strains  and  pressures  relevant  to  the 
laboratory  membrane,  £  is  well-approximated  by  the  truncated  series 


£  =  dip  +  dzp 3  +  ds5p5,  dss  =  1  —  d\  —  d3, 


(7.9) 
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while  p  can  be  approximated  by 

P  =  h  i  +  4-  /ss^5,  fss  =  1  fi  /3-  (7.10) 

Once  the  series  (7.7)  for  p  is  determined  to  the  desired  accuracy,  the  result  can  be  substituted  for  p  —  F -1(f)  in 
the  power  series  representation  (6.21)  of  W(p ),  yielding  according  to  (7.5)  a  series  solution  for  the  equation  of  the 
true  middle  surface: 

w(r)  =  a  {  u50  +  C  c2n  [F~l{r/a)]2n  },  (7-H) 

n= i 

or,  in  dimensionless  form: 

w(t)  =  Wo  +  C  C2n  (7-12) 

n  =  1 

and  we  note  that  since  F~1(  1)  =  1,  we  have  u5(l)  =  0,  as  required.  In  Figure  4  we  show  a  rather  extreme  case  (zero 
pre-strain  and  an  /-number  of  2)  of  the  difference  between  the  true  surface  and  the  surface  of  lateral  displacements. 
For  pre-strains  of  as  little  as  0.1%  and  high  /-numbers,  say  // 8,  the  difference  is,  however,  negligibly  small  (on  the 
order  of  tenths  of  a  wavelength) . 


Deviation  of  True  Surface 
from  Lateral  Displacement  Surface 


-1.0  -0.5  0.0  0.5  1.0 

?  =  r/a 


FIG.  4.  Deviation  of  True  Surface  from  Surface  of  Lateral  Displacements  (// 2  Membrane) 
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VIII.  THE  MEMBRANE  SURFACE  AS  A  FUNCTION  OF  F-NUMBER 


Referring  to  Figure  3,  we  note  that  for  r  «  a,  i.e.,  near  the  vertex  of  the  membrane  where  r  =  0,  the  equation  of 
the  surface  is  approximated  by  the  first  two  terms  of  (7.11),  viz., 


u(r)  «  a  l^o  +  Ccifl  (^)'i 


C/i2  2 

Wo - r% 

a 


(8.1) 


where  we  replaced  C2  =  — 72  =  —  1  in  the  final  equality.  This  is  recognized  to  be  the  equation  of  a  paraboloid  with 
origin  at  O ,  having  a  focal  length  /  given  by 


/  = 


a  _  a  d\ 

ACfl  ~  Tc' 


(8.2) 


where  the  final  equality  follows  from  the  first  relation  in  (7.8).  In  all  that  follows,  we  consider  this  “paraxial”  focal 
length  to  be  the  focal  length  of  the  membrane.  It  is  convenient  to  replace  /  by  the  membrane  /-number ,  //#,  defined 

by 


//#  = 


2a1 


where  2a  is  the  membrane  diameter.  From  (8.2)  and  (8.3)  it  then  follows  that 


//# 


8  C' 


(8.3) 


(8.4) 


which  relates  the  /-number  to  the  mechanical  properties  of  the  membrane  (recall  that  d\  depends  on  both  e  and  v 
according  to  (7.6),  while  C  depends  on  q  —  pa/Eh  according  to  (6.24)  and  (4.22)  ). 

In  order  to  determine  the  equation  of  the  membrane  surface  for  a  given  /-number  //#,  and  given  values  of  e  and 
v,  we  begin  by  noting  from  (8.4)  that  C  must  in  this  case  be  given  by 


_  <*i  _  [(1  ~  0  +  ^0  ( 1  ~  V)Y 

8//#  "  8//# 


(8.5) 


where  the  last  expression  follows  from  (7.6).  From  the  definition  (6.24)  we  obtain  the  following  relation  between  C , 
60,  and  the  variable  q: 


q  —  4boC. 


Substitution  of  this  expression  into  (6.13)  yields 


B  = 


2\/  6q 


(8.6) 


(8.7) 


(8.8) 


These  expressions  for  C  and  B  as  functions  of  bQ  are  now  used  in  (6.17),  to  determine  60  as  a  solution  of 

[(1  -  v)  +  (3  -  v )62  +  (5  -  v)bA  +  (7  -  v)b6  +  (9  -  v)b8  +  ...]  -  ^  =  °, 

for  a  given  value  of  the  pre-strain  e.  Once  bo  is  determined,  the  coefficients  in  the  power  series  for  the  displacement 
components  and  other  variables  can  be  computed.  The  pressure  for  this  value  of  bo-  and  the  given  values  of  the 
parameters  e,  v,  and  //#,  is  found  from  the  definition  (4.22): 


P  =  q 


Eh 


(8.9) 


after  q  has  been  computed  from  (8.6).  The  maximum  displacement  is 

wo  =  a  wq, 


(8.10) 


where  w0  has  been  computed  according  to  (6.27). 
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IX.  ESTIMATION  OF  THE  PRE-STRAIN  AS  A  FUNCTION  OF  F-NUMBER 


Suppose  that  at  a  given  pressure  p  the  membrane  is  characterized  as  having  a  particular  /-number  //#.  The 
amount  of  pre-strain  e  for  these  values  of  p  and  //#  is  difficult  to  measure  directly,  but  can  be  calculated  as  follows. 
From  (8.5),  (8.6),  and  (4.22),  we  have 


V  = 


[1  +  60(1  -  v) 
2//# 


(9.1) 


which  can  be  solved  for  e  as  a  function  of  the  known  parameters,  and  the  unknown  coefficient  fco'- 


■  =  ‘  +  (9'2) 

Since  p  is  known,  we  can  compute  q  from  (4.22),  and  then  from  (6.24)  and  (6.13)  get  C  and  B  as  functions  of  the 
unknown  coefficient  6q: 


C  = 


9 

460  ’ 


(9.3) 


These  expressions  for  e  and  B  as  functions  of  bo  are  now  substituted  in  (8.8),  which  is  solved  for  bo.  This  value  of  bo 
can  then  be  used  in  (9.2)  to  obtain  the  desired  estimate  of  the  pre-strain  €  at  the  given  values  of  p  and  //#. 

We  also  note  that  for  e/[l  +  60  ( 1  -  v )]  <<  1,  we  obtain  from  (9.1)  a  linear  approximation  of  p  as  a  function  of  e: 

P  »  ——  7T7T77  [  1  +  Ml  -  v)]2  -  —7 tM1  +  M1  -  v)\*  =  Po  -  me,  (9.4) 

a  2//#  a  //# 

where  the  p-intercept  and  slope  m  are  given  by 

Po  =  IT  277#^  +  “  I/)]2’  and  m  =  +  b°(l  ~  ^9'5^ 

The  approximation  (9.4)  yields  the  following  linear  approximation  for  the  pre-strain  as  a  function  of  pressure,  for  a 
given  //#: 


m 


(9.6) 


X.  ESTIMATION  OF  THE  PRE-STRAIN  AS  A  FUNCTION  OF  MAXIMUM  LATERAL 

DISPLACEMENT 


It  has  been  found  useful  to  have  at  least  an  approximate  solution  of  the  following  problem:  what  pre-strain  e  is 
required  to  obtain  a  desired  maximum  lateral  displacement  of  the  membrane  for  a  given  pressure  pi  We  note  that 
from  (6.21)  and  (6.27)  we  have  the  following  expression  for  the  maximum  displacement  Wo: 

W(0)  =  Wo  =  a  wo  —  —qC  ^2  c2n, 

n=l 


or,  using  (6.25): 

Wq  =  aC  (72  +  74  B2  +7 qBa  +  78  B6  +  . . .  H-  720  #18 )» 

where  we  recall  from  (6.26)  that  7 2n  is  the  numerical  coefficient  of  B2n~2  in  each  series  coefficient  C2n>  e.g.,  72  =  1, 
74  =  1/2,  76  =  5/9,  and  7s  =  55/72.  Using  (6.24)  to  replace  C  yields 


The  coefficient  B  depends,  in  general,  on  bo  and  q  according  to  its  definition  (6.13),  and  we  recall  that  bo  is  determined 
for  a  given  q,  v,  and  e  by  solving  (6.17). 

However,  for  e  =  0,  equation  (6.17)  reduces  to 

1  -  v)  +  (3  -  v)bo  +  (5  -  v)bA  +  (7  -  v)h  +  (9  -  u)ba  +  ...  =  0.  (10.2) 

which  can  be  solved  for  B  without  regard  for  q,  requiring  only  a  given  value  for  v  [the  coefficients  62„,  n  >  1,  depend 
only  on  B,  according  to  (6.15)].  Denoting  by  B0  the  value  of  B  corresponding  to  e  =  0,  we  have  from  (6.13): 

*»  - 1  (0'3-  <“-3> 

for  the  value  of  60  when  e  =  0.  This  can  be  substituted  for  b0  in  (10.1)  to  obtain  the  maximum  displacement  with 
zero  pre-strain: 

Wo  =  aqx^B^3  (72  +  74  +  76  B4  +  7s  Be  +  - . .  +  720  #18)  =  JCaq1^,  (lO-4) 

where  we  have  introduced  a  new  constant  /C  defined  by 


K  =  Bl/3  (72  +  74  B2  +  76  B4  +  7s  B6  +  . . .  +  720fl18  )•  (10-5) 

We  remark  that  our  equation  (10.4)  is  equivalent  to  Campbell’s  equation  (25),  where  JC  —  0.653  corresponds  to  his 
use  of  a  Poisson’s  ratio  v  =  0.3.  We  have  verified  this  value  of  K  with  our  computations,  and  find  that  when  we 
increase  v  to  v  -  0.4,  the  value  decreases  to  K  -  0.626.  Once  this  constant  has  been  determined  for  a  given  value  of 
i/,  we  can  solve  (10.4)  for  the  pressure: 


Po 


__  Eh  /W0y 
a  V  JCa) 


(10.6) 


which  is  the  pressure  required  to  produce  a  desired  maximum  displacement  Wo  when  there  is  zero  pre-strain. 

Now,  in  the  opposite  limit  of  large  pre-strain,  we  find  numerically  that  the  parameter  B  tends  to  a  number  much 
less  than  1,  hence  the  terms  in  (10.1)  involving  B  can  be  neglected  to  yield  the  approximation 


Wo  =  %  jh,  (10-7) 

4  60oo 

where  we  denote  by  b0oa  the  value  of  b0  for  large  values  of  e.  In  the  same  approximation,  since  the  coefficients  62„, 
n  >  1,  are  proportional  to  powers  of  B,  the  determining  equation  (6.17)  reduces  to  simply 


1  —  v  —  — =  0  =>  booo  =  - - •  (10.8) 

bocc  1  -  " 

Substituting  (4.22),  and  this  expression  for  bQca,  in  (10.7)  yields  the  following  expression  for  the  pressure  under  large 
pre-strain  conditions: 


Poo  — 


Wo 

- e. 

1  -  v 


(10.9) 


The  curve  of  pressure  versus  pre- strain  for  a  given  value  of  Wo  has  been  found  to  be  well- approximated  by  the 
linear  function  obtained  from  the  sum  of  the  two  extreme  cases  (10.6)  and  (10.9): 


p(c)  «  po  +  rne, 


(10.10) 


where  the  p-intercept  pQ ,  and  slope  m,  are  given  by 


Po 


Eh  /Wo\3 
a  \Ka) 


and 


m  = 


,  Eh  Wo 
4——v' 


(10.11) 


respectively.  Thus,  by  measuring  the  pressure  p  required  to  produce  a  (measured)  maximum  displacement  Wo,  we 
can  invert  (10.10)  to  obtain  an  estimate  for  the  amount  of  initial  strain  existing  in  the  membrane: 


p  -  Po 

€  «  - . 

m 

Unfortunately,  a  precise  direct  measurement  of  Wo  is  also  difficult  to  obtain. 


(10.12) 
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XI.  COMPARISONS  WITH  REFERENCE  SURFACES 


In  Figure  5  we  have  illustrated  again  the  displacement  of  point  P  of  an  initially  flat  membrane  to  point  Pf  of  its 
equilibrium  deformed  configuration.  We  have  also  shown  there  a  sphere  of  radius  Rc  with  center  at  C ,  which  has  been 
chosen  to  coincide  with  the  membrane  at  the  vertex  V  and  along  the  circular  membrane  edge  r  —  a.  The  equation  of 
the  middle  surface  of  the  deformed  membrane,  in  a  frame  with  origin  at  O,  is  given  by  (7.5),  i.e., 

z  =  tn(r).  (11.1) 

The  center  0  of  the  middle  plane  Z  -  0  is  displaced  under  the  deformation  to  the  vertex  V  of  the  membrane,  which 
has  coordinates  (0,  u;0),  where  wq  =  w( 0)  =  W{ 0)  =  Wo  is  the  maximum  lateral  displacement. 


In  Figure  5,  a  point  with  Cartesian  coordinates  {z,  y ,  z }  in  the  frame  with  origin  at  O  has  coordinates  {a?,  y,  z}  in 
a  frame  with  origin  at  the  vertex  V.  The  two  axial  coordinates  are  related  by  a  translation  through  the  maximum 
displacement  wo,  i.e., 

z  ~  z  —  wo,  (11 .2) 


which  assures  that  the  vertex  V  defined  by  z  =  0  corresponds  to  z  =  w0  with  respect  to  O.  The  equation  of  the 
membrane  surface  in  a  frame  with  origin  at  V  is  thus 

z  =  w(r)  —  wo.  (11.3) 


In  this  coordinate  system,  the  equation  of  a  paraboloid  with  focal  length  fp  has  the  particularly  simple  form 


while  points  on  the  surface  of  the  upper  hemisphere  of  a  sphere  of  arbitrary  radius  Rc  satisfy  the  equation 

-  -Rc  +  V&c  "  r2, 


(11.4) 


(11.5) 


which  is  easily  derived  from  the  geometry  of  Figure  5.  We  note  that  for  r  <<  Rc  (  i.e.,  near  the  vertex  V  ),  the 
equation  of  the  sphere  is  approximated  by 
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(11.6) 


which  we  obtained  by  keeping  only  the  first  two  terms  in  the  binomial  expansion  of  the  square  root  function  in  (11.5). 

According  to  (11.4),  this  represents  a  paraboloid  of  focal  length  fs  =  Rc/ 2.  Translating  back  to  the  coordinate 

system  with  origin  at  the  membrane  center  O,  the  equations  for  spherical  and  paraboloidal  surfaces  passing  through 
the  vertex  V  are  given  by 

zs  —  ws{r)  -  w0  -  Rc  +  yj R2C  -  r2  ,  (H-7) 

«  Wo  -  =jr-r2,  r  «  Rc,  (11.8) 

Zrtc 


and 


zp  =  wp(r)  =  wo  -  4j^r2> 


(11.9) 


where  the  approximation  (11.8)  holds  near  the  vertex. 

We  note  that  the  plane  z  —  0  intersects  the  membrane  surface  in  a  circle  of  radius  a,  which  suggests  choosing  as 
a  reference  sphere  one  obtained  by  requiring  that  the  plane  2=0  also  cut  the  spherical  surface  in  a  circle  of  radius 
r  =  a,  as  shown  in  Figure  5.  From  (11.7)  it  follows  that  the  maximum  displacement  w o,  the  membrane  radius  a,  and 
the  radius  Rc  of  the  sphere,  must  then  satisfy  the  condition 


wo  =  Rc  -  \I  R%  —  a2  • 


(11.10) 


This  can  be  solved  for  the  radius  of  curvature,  to  obtain 


Rc 


a2  +  Wp 

2w0 


(11.11) 


A  reference  paraboloid  can  be  defined  by  the  same  condition,  viz.,  that  the  plane  z  —  0  cut  its  surface  in  a  circle  of 
radius  r  —  a.  For  this  to  occur,  we  must  have  from  (11.9): 


a2  ,  q2 

w0  =  77-  =>  fp  =  — • 

4/p  4iuo 


(11.12) 


From  (11.11)  and  (11.12),  it  follows  that  for  reference  surfaces  satisfying  these  criteria,  the  radius  of  curvature  of  the 
sphere  and  the  focal  length  of  the  paraboloid  must  be  related  by  the  expression 


Rc  =  2/p  +  y-  =  2(/p  +  ^)  =  2/s,  where  fs  =  fp  +  (11.13) 


Thus,  in  the  approximation  r  «  Rc,  this  reference  sphere  is  approximated  not  by  a  paraboloid  of  focal  length  fp, 
but  rather  by  one  of  longer  focal  length  fs  ■  Both  focal  lengths  are  determined  by  the  maximum  displacement  w0 
which,  for  a  given  membrane  /-number,  is  obtained  from  the  algorithm  leading  to  equation  (8.10).  It  should  be  noted 
that  in  general  neither  of  these  focal  lengths  is  the  same  as  the  membrane  focal  length  /  defined  by  (8.2). 

The  last  observation  suggests  a  different  set  of  reference  surfaces,  i.e.,  we  can  choose  a  reference  paraboloid  having 
focal  length  fp  =  /,  and  a  reference  sphere  having  a  radius  of  curvature  that  is  twice  this  focal  length,  i.e.,  Rc  —  2f. 
With  these  choices,  regions  near  the  vertex  of  any  one  of  the  three  surfaces  will  have  the  same  focal  length. 

For  either  set  of  reference  surfaces,  we  can  compute  the  deviation  Ato(r)  of  the  membrane  surface  from  the  spherical 
and  paraboloidal  reference  surfaces,  i.e., 


A ws(r)  =  w(r)  -  ws{r)  =  w(r)  -  w0  +  Rc  -  y/ R%  —  r2 


(11.14) 


and 

(11.15) 

respectively.  Since  these  differences  are  typically  on  the  order  of  microns,  it  is  convenient  to  graph  them  in  units  of 
the  wavelength  A  of  light  used  to  probe  the  membrane  mirror,  i.e.,  we  plot 


Awp(r)  =  w(r)  —  wp(r )  —  w(r) 


w0  +  77 ~r 
AfP 
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(11.16) 


and  An,„(r,A)  = 

As  an  example,  Figure  6  is  a  graph  of  — Au>s(£,A)  as  a  function  of  £  for  an  f/S  membrane  with  three  different 
levels  of  pre-strain,  where  we  have  used  the  reference  sphere  that  contacts  the  membrane  at  both  the  vertex  and  the 
membrane  edge.  For  0.25%  pre-strain,  the  graph  indicates  a  nearly  95%  reduction  in  the  maximum  deviation  from 
a  sphere,  which  would  imply  a  similar  reduction  in  the  spherical  aberration  of  the  wavefront  reflected  from  an  f/S 
membrane  mirror. 


- P*rc*nt  Pr#-S  train  •  0.0 

— —  P«rc#nt  Pr«-8  train  •  0.1 


FIG.  6.  Deviation  of  f/S  Membrane  from  Reference  Sphere 


We  are  usually  interested  in  this  maximum  deviation  from  a  reference  surface,  which  occurs  at  a  value  of  r  for 
which 

-^-[Au>s(r)]  =0  or  —  [Atup(r)]  =  0.  (11.17) 

dr  dr 

To  determine  the  left-hand  sides  of  these  expressions,  we  first  use  the  chain  rule  to  compute  dw/dr ,  i.e.,  we  have 

r 

w(r)  —  aiD(£(r)),  where  £(r)  = 


hence 


However, 


dw  _  dw  d£  _  dw 
dr  a  d£  dr  d£ 

w(0  =  W(p(£)), 


so  the  chain  rule  can  be  applied  again  to  obtain 


dw  __  dW  dp 
dr  dp  d£ 


Recalling  from  (7.7)  that 


P  =  hi  +  h?  +  h?  +  ■■■  , 


(11.18) 


(11.19) 
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(11.20) 


we  substitute  this  result  together  with  (6.22)  into  (11.18),  yielding 

^  =  2  ( /i  +  3  M2  +  5/s£4  +  ...)- £  nc*nP2n, 
dr  n  — 


p  s 


where  we  made  the  replacement  u>2n  =  C  C2n  •  It  is  easy  to  show  that 


dws 

dr 


\pk  -  & 


■ — ■  _RC 

where  #c  =  — , 
a 


and 


dwp 

dr 


2/p 


where  fp  =  —, 
a 


hence  the  value  of  r  or  £  =  r/a  that  will  give  the  maximum  deviations  must  be  a  solution  of  either 


2/Wl  -  i2  (fi  +  3 M2  +  5/s£4  +  ...)C£nc2„p2"  +  />£  =  0, 


(11.21) 


(11.22) 


(11.23) 


n=l 


4  fp  ( /1  +  3  fs^2  4-  5 /s  £4  +  . . .)  C  n C2n  p2n  +  p£  —  0-  (11.24) 

n=l 

A  graph  of  the  maximum  deviation  of  //4,  //6,  and  // 8  membranes  from  a  reference  sphere  as  a  function  of  percent 
pre-strain  is  shown  in  Figure  7,  below.  The  marked  improvement  (over  the  zero  pre-strain  case)  to  approximately  six 
waves  of  deviation  requires  only  0.25%  pre-strain  for  an  // 8  membrane,  but  more  than  2.0%  pre-strain  for  an  // 4 
membrane. 


FIG.  7.  Maximum  Deviation  from  Reference  Sphere- versus-Pre- Strain  (// 4,  // 6,  and  // 8  Membranes) 
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XII.  RAY  DEVIATIONS  FROM  A  MEMBRANE  MIRROR  ILLUMINATED  BY  AN  ON-AXIS  POINT 

SOURCE 

The  solution  (11.1)  for  the  deformed  membrane  surface  defines  a  surface  of  revolution.  The  equation  of  such  a 
surface  can  also  be  written  as 


<£(r,z)  =  0. 


(12.1) 


where  the  function  <j>  is  defined  by 


<^(r,  z)  =  w(r)  —  z. 


(12.2) 


The  value  of  the  gradient  of  this  function  at  a  point  of  the  membrane  with  coordinates  (r,  z ),  i.e.,  V<£(r,  z)}  is  a  vector 
normal  to  the  membrane  surface  at  that  point,  hence  the  unit  normal  vector  at  any  such  point  is  given  by 


V* 

"  I v*r 


(12.3) 


Thus,  from  (12.3)  and  (12.2),  we  have  for  the  unit  normal  to  a  surface  of  revolution,  expressed  in  terms  of  the 
orthonormal  basis  {er ,  e$ ,  ez } : 

n(r)  =  —  — — ,  where  V(r)  =  \J  [w'(r)]2  +  1  ,  (12.4) 


and  the  prime  denotes  a  derivative  with  respect  to  r. 

Suppose,  now,  that  an  incident  ray  propagating  in  a  direction  specified  by  unit  vector  fe*  strikes  a  point  P  of  the 
surface  as  illustrated  in  Figure  8.  At  this  point  of  incidence  there  will  originate,  in  general,  both  a  reflected  and  a 
transmitted  (or  refracted)  ray.  The  transmitted  ray  is  of  no  interest  here  (for  a  high  reflectance  surface,  the  transmitted 
ray  is  essentially  nonexistent).  From  Figure  3  and  the  (geometrical)  definitions  of  the  dot  and  cross  products  of  two 
vectors,  we  have 


n-fef  =  cos  (tt  —  &i)  =  —  cos  0i,  (12.5) 

nxki  ~  sin(7r  —  9i)e$  =  sin  e$ ,  (12.6) 


where  0,-  is  the  angle  of  incidence. 


FIG.  8.  Geometry  of  Incident  and  Reflected  Rays 
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The  direction  of  the  reflected  ray  is  specified  by  a  unit  vector  kr.  To  facilitate  the  computation  of  this  vector  we 
introduce  at  P  a  unit  vector  v  tangent  to  the  meridional  curves  on  the  surface,  defined  by 

v  =  nXee-  (12.7) 

The  unit  ray  vectors  ki  and  fer  lie  in  a  plane  called  the  plane  of  incidence,  spanned  by  the  orthonormal  pair  {n,v}. 
Thus,  the  unit  ray  vectors  can  be  written  in  terms  of  components  along  n  and  v  as 

ki  =  —cos din  -  sinOiV,  (12.8) 

kr  =  cos02-n  -  sin#;*;,  (12.9) 

where  we  used  the  law  of  reflection,  which  states  that  the  angle  of  reflection  equals  the  angle  of  incidence  0*.  By 
subtracting  (12.8)  from  (12.9)  and  using  (12.5)  in  the  result,  we  obtain  fcr  -  ki  =  2  cos  0*71  =  -2 (n-fe,)n,  from 
which 


kr  =  ki  -  2 (n*fe*)n.  (12.10) 

Equation  (12.10)  is  a  coordinate-free  expression  for  the  direction  of  the  reflected  ray  in  terms  of  the  given  incident 
direction  and  the  calculated  surface  unit  normal  vector. 

Referring  now  to  Figure  9,  we  consider  rays  emanating  from  a  point  source  located  at  C,  the  center  of  curvature  of 
the  membrane  (not  the  reference  sphere,  as  was  the  case  in  the  previous  Section),  which  has  position  vector 

L  =  -Le2,  1  =  2/ -wo,  (12.11) 

with  respect  to  O,  as  can  be  easily  deduced  from  the  Figure.  Any  such  ray  propagates  radially  outward  from  C,  and 
we  show  a  particular  ray  striking  the  surface  at  an  arbitrary  point  P,  which  is  then  reflected  to  some  point  Q. 


We  denote  by 

x  —  reT  4-  zez  —  r  er  +  w(r)ezi  where  z  =  u>(r),  (12.12) 

the  position  vector  of  P  with  respect  to  O,  and  by  xj  the  position  vector  of  Q  with  respect  to  O,  so  that 

Xd  —  x  4-  dkr  —  x  4-  d[ki  —  2  (n*fe2-)n],  (12.13) 

where  d  —  \xd  —  ®|  is  the  distance  between  P  and  Q  along  the  line  parallel  to  fer,  and  (12.10)  was  used  to  replace 
kr  in  the  final  equality. 
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The  position  vector  of  Q  with  respect  to  C  is  called  the  ray  deviation  vector ,  which  we  denote  by  A.  From  Figure 
9  we  have 

A  =  xd-  L,  (12.14) 

hence  from  (12.13)  and  (12.14)  it  follows  that 

A  =  x  -  L  +  d[ki  -  2  (n-ki)n}.  (12.15) 

For  convenience,  we  introduce  a  new  variable  £  defined  by 

£  =  z  +  L  =  it;(r)  +  L,  (12.16) 

so  that 

x-L  =  rer  +  zez+Lez  =  rer+£ez.  (12.17) 

Now,  the  incident  ray  unit  vector  ki  in  this  case  is  just  a  unit  vector  along  the  position  vector  of  P  with  respect  to 
C,  i.e., 

ki  =  - - ~~r  =  ^  — ,  where  Pc(r,£)  =  \/ r2  4-  £2  •  (12.18) 

\x  —  L  |  Pc 

Substituting  (12.18)  in  (12.15)  yields 


A  =  ®*L  +  —  {  (a?  -  L)  -  2  [n-(  x  -  L)]n }. 
Pc 

The  radial  and  lateral  components  of  the  ray  deviation  vector  are  thus 


(12.19) 


4  A  d  ,  2wf  .  /  r  v , 

Ar  =  er-A  =  r  +  —  {r  -  —  [»•(*  “  £)], 


Az  =  ez-A  =  C  +  ^{C  +  ^ [«•*(*  -  •£)]> 
or,  since  w(x  —  T )  —  (ru/  -  £  )/P,  we  have 

Ar  =  r  +  A[r  -  ^■(ro'  -  C)L  A2  =  C  +  ^  [C  +  ( rw>  ~  01- 

Factoring  1/P2  from  the  square-bracketed  terms,  recalling  that  P2  =  ( w ')2  +  1,  these  can  also  be  written  as 

Ar  =  r  +  2v^{Cw>  +  m)’  Az  =  c  +  2v^{rw'  ~  Ca)’ 

where  we  have  introduced  for  convenience: 

O'  =  i[l  -  K)2]. 

The  factor  d/V^Vc  can  be  eliminated  between  the  two  equations  of  (12.20),  yielding  the  following  relation  between 
the  ray  deviation  components: 


(12.20) 


(12.21) 


_  £ w'  +  ra  (r2  -  £2)^  -  2r£a 

r  r  w'  —  £  Q'  z  r  w'  —  £  Q 


(12.22) 


For  a  given  value  of  the  longitudinal  deviation  Az ,  we  can  use  (12.22)  to  graph  the  transverse  deviation  Ar  as  a 
function  of  the  coordinate  r  of  the  point  P  on  the  membrane  from  which  the  ray  was  reflected. 

It  is  also  useful  to  compute  the  angular  deviation  8  of  a  reflected  ray  from  the  incident  ray.  This  deviation  is  just 
twice  the  angle  of  incidence,  i.e.,  8  =  2  9{.  From  (12.5),  (12.6),  (12.4),  and  (12.18),  we  have 

ri'ki  —  —  cos  Oi  —  (  r  w'  —  £),  (12.23) 


PPc 

nxki  =  sin  Oi  eg  —  {(w*  +  r)e$, 

PPc 


hence  the  ray  deviation  angle  is  given  by 

4  =  Sft=2i«-(£2^). 

\rw'  —  C  / 


(12.24) 

(12.25) 
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XIII.  COMPARISONS  WITH  LABORATORY  MEASUREMENTS 


Two  types  of  measurements  were  made  [10]  on  an  11-inch  28-cm)  diameter  laboratory  membrane  of  125//  thick 
Upilex,  having  a  Young’s  modulus  of  E  =  8800  N/m 2  and  a  Poisson’s  ratio  of  v  =  0.4.  In  the  first  type  of  measurement, 
the  annulus  vacuum  pressure  was  fixed  at  some  value  (po  =  0  inches  of  water  in  one  case,  and  po  =  20  inches  of  water 
in  the  other  case),  and  the  central  displacement  (sag)  wq  was  measured  as  the  inner  membrane  vacuum  pressure  p 
was  increased.  A  measurement  of  inner  membrane  pressure  was  recorded  at  every  0.010-inch  (10  mil)  increment  of 
the  central  sag  as  shown  in  the  following  tables: 


po  =  0  inches  water 


wo  (mm) 

p  (inches  water) 

0.254 

0.21 

0.508 

0.42 

0.762 

0.63 

1.016 

0.85 

1.270 

1.09 

1.524 

1.30 

1.778 

1.58 

2.032 

1.87 

2.286 

2.19 

2.540 

2.50 

2.794 

2.83 

3.048 

3.25 

3.302 

3.64 

3.556 

4.12 

3.810 

4.60 

4.064 

5.13 

4.318 

5.73 

4.572 

6.30 

4.826 

7.03 

5.080 

7.77 

po  =  20  inches  water 


wq  (mm) 

p  (inches  water) 

0.254 

0.32 

0.508 

0.67 

1.016 

1.40 

1.524 

2.07 

2.032 

2.85 

2.54 

3.70 

3.048 

4.54 

3.556 

5.65 

4.064 

6.78 

4.572 

8.08 

5.08 

9.62 

For  each  data  point  (iu0,p),  the  Hencky-Campbell  theory  was  used  to  compute  the  predicted  pre-strain  e  existing  in 
the  inner  membrane.  These  computed  values  are  plotted  versus  the  measured  central  sag  in  the  following  two  graphs. 
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Computed  Pre-Strain  -vs-  Measured  Central  Sag 
(Annulus  Pressure  =*  20"  water) 
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These  graphs  show  that  although  the  annulus  pressure  was  maintained  at  a  nearly  constant  value,  the  amount  of 
strain  in  the  inner  membrane  attributable  to  the  “pre-strain”  occurring  in  the  Hencky- Campbell  theory  decreases 
as  the  inner  membrane  vacuum  pressure  is  increased.  This  is  not  surprising  since  the  laboratory  membrane  is  not 
clamped  on  the  inner  ring,  hence  an  increase  in  inner  membrane  pressure  can  cause  some  of  the  membrane  material 
over  the  annulus  chamber  to  slip  across  the  inner  ring  to  become  part  of  the  inner  membrane.  The  effect  of  this 
membrane  migration  is  a  reduction  of  the  total  strain  in  the  inner  membrane  (compared  to  the  strain  that  would  exist 
if  the  membrane  had  been  clamped  on  the  inner  ring),  which  manifests  in  the  model  as  a  reduction  of  the  effective 
pre-strain.  The  next  two  graphs  are  plots  of  the  inner  membrane  pressure  against  central  sag,  where  the  points  are 
the  data  points  from  the  tables  on  the  last  page,  and  the  curves  reflect  model  predictions  using  the  mean  value  (not 
a  particularly  good  estimator)  of  the  computed  pre-strains  shown  in  the  previous  two  graphs. 


-  Model  estimate  using 
average  pre-strain 
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•  i  oh  data 

- Mode!  estimate  using 


In  the  second  type  of  measurement  the  inner  membrane  vacuum  pressure  was  adjusted  until  the  central  deflection 
was  wo  =  2.1  mm.  This  was  done  for  two  different  annulus  vacuum  pressures.  The  inner  membrane  vacuum  pressure 
required  for  the  low  annulus  pressure  case  was  p  =  2.0  inches  of  water,  while  the  high  annulus  pressure  case  required 
a  pressure  of  p  =  8.1  inches  of  water.  These  data  were  used  to  calculate  pre-strains  of  0.05%  and  0.24%,  respectively. 
The  predicted  deviations  of  the  membrane  surfaces  from  a  reference  sphere  (contacting  the  membrane  at  its  vertex 
and  along  the  inner  ring)  in  the  two  examples  is  shown  in  the  next  graph.  The  effect  of  a  higher  pre-strain  on  the 
shape  is  rather  dramatic,  causing  the  maximum  deviation  to  be  reduced  from  over  22  waves  to  about  6  waves. 


- P  .  2.0"  water  (f/#  «  8.55,  e  -  0.05%) 

- P  -  8.1  ”  water  (f/*  =  8.39,  e  -  0.24%) 


In  each  case  the  membrane  mirror  was  being  illuminated  by  laser  light  from  a  pinhole  (simulated  point  source)  located 
at  the  center  of  curvature  of  the  membrane.  Measurements  were  made  of  the  minimum  geometric  spot  size  of  the 
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light  reflected  from  the  mirror.  The  measured  minimum  spot  size  was  roughly  7-8  mm  for  the  low  pre-strain  case, 
and  3-4  mm  for  the  higher  pre-strain.  In  the  next  graph  we  show  plots  of  the  computed  transverse  ray  deviations  in 
the  z-planes  where  the  maximum  transverse  ray  deviations  were  minimized  (z  —  35.7  cm  and  z  =  9.2  cm  from  the 
membrane  center  of  curvature  for  the  low  and  high  pre-strain,  respectively). 


Computed  Ray  Deviation  -vs-  Normalized  Aperture  Radial  Coordinate 
(Fixed  2.1  mm  Centra!  Sag) 

- P-  2.0"  water,  Z  -35.7  cm  from  Center  of  Curvature 


The  predicted  peak-to-valley  ray  deviations,  which  correspond  directly  to  the  spot  diameters,  are  seen  to  be  on  the 
order  of  7  mm  and  2  mm,  respectively,  somewhat  smaller  than  the  measured  spot  diameters.  This  is  to  be  expected, 
since  our  axially  symmetric  model  does  not  account  for  asymmetric  aberrations  observed  in  the  laboratory  membrane 
(especially  astigmatism),  which  contribute  to  the  larger  measured  spot  sizes. 


[1]  Campbell,  J.  D.,  Quart.  J.  Mech.  and  Appl.  Math.  IX,  84-93,  1956. 

[2]  Fichter,  W.  B.,  NASA  Technical  Paper  3658,  July  1997. 

[3]  Fung,  Y.  C.,  Foundations  of  Solid  Mechanics,  Prentice-Hall,  Englewood  Cliffs,  1965. 

[4]  Greschik,  G.,  Mikulas,  M.  M.,  and  Palisoc,  A.,  Proceedings  of  the  39th  Structures,  Structural  Dynamics,  and  Materials 
Conference  and  Adaptive  Structures  Forum,  Long  Beach,  CA,  April  20-23,  1998.  Volume  4,  2772-2782,  AIAA.  AIAA-98- 
2101-CP. 

[5]  Hendcy,  H.,  Zeitschrift  fiir  Mathematik  und  Physik  63,  311-317,  1915. 

[6]  Marker,  D.  K.,  Carreras,  R.  A.,  and  Jenkins,  C.  H.,  Second  International  Symposium  on  Intelligent  Automation  and 
Control,  World  Automation  Congress  (WAC),  9-14  May  1998,  Anchorage,  Alaska. 

[7]  Marker,  D.  K.,  Carreras,  R.  A.,  Wilkes,  J.  M.,  and  Jenkins,  C.  H.,  Ninth  Annual  Conference  Laser  Optics  98,  Vavilov 
State  Optical  Institute,  22-26  June  1998,  St.  Petersburg}),  Russia. 

[8]  Press,  W.  H.,  Flannery,  B.  P.,  Teukolsky,  S.  A.,  and  Vetterling,  W.  T.,  Numerical  Recipes:  The  Art  of  Scientific  Computing, 
Cambridge  University  Press,  New  York,  1986. 

[9]  Rotge,  J.  R.,  Marker,  D.  K.,  Carreras,  R.  A.,  Jenkins,  C.  H.,  Wilkes,  J.  M.,  Duneman,  D.,  SPIE,  20-27  July  1998,  San 
Diego,  California. 

[10]  Rotge,  J.  R.,  private  communication. 

[11]  von  Karman,  T.,  Encyklopddie  der  Mathematischen  Wissenschaften  IV,  349,  1910. 


28 


DISTRIBUTION  LIST 


AUL/LSE 

Bldg  1405  -  600  Chennault  Circle 

Maxwell  AFB,  AL  361 12-6424  1  cy 

DTIC/OCP 

8725  John  J.  Kingman  Rd,  Suite  0944 
Ft  Belvoir,  VA  22060-6218  2  cys 

AFSAA/SAI 

1580  Air  Force  Pentagon 

Washington,  DC  20330-1580  1  cy 

AFRL/PSTL 

Kirtland  AFB,  NM  87 1 1 7-5776  2  cys 

AFRL/PSTP 

Kirtland  AFB,  NM  871 17-5776  1  cy 

AFRL/DE/Dr.  Good 

Kirtland  AFB,  NM  87117-5776  1  cy 

Official  Record  Copy 

AFRL/DEB S/J ames  M.  Wilkes  1 0  cys 


29 


